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This planetary boundaries framework update finds that six of the nine boundaries are transgressed, suggesting poeaes 


that Earth is now well outside of the safe operating space for humanity. Ocean acidification is close to being 
breached, while aerosol loading regionally exceeds the boundary. Stratospheric ozone levels have slightly re- 
covered. The transgression level has increased for all boundaries earlier identified as overstepped. As primary 
production drives Earth system biosphere functions, human appropriation of net primary production is pro- 
posed as a control variable for functional biosphere integrity. This boundary is also transgressed. Earth 
system modeling of different levels of the transgression of the climate and land system change boundaries il- 


lustrates that these anthropogenic impacts on Earth system must be considered in a systemic context. 


INTRODUCTION 
The planetary boundaries framework (1, 2) draws upon Earth 
system science (3). It identifies nine processes that are critical for 
maintaining the stability and resilience of Earth system as a 
whole. All are presently heavily perturbed by human activities. 
The framework aims to delineate and quantify levels of anthropo- 
genic perturbation that, ifrespected, would allow Earth to remain in 
a “Holocene-like” interglacial state. In such a state, global environ- 
mental functions and life-support systems remain similar to those 
experienced over the past ~10,000 years rather than changing into a 
state without analog in human history. This Holocene period, which 
began with the end of the last ice age and during which agriculture 
and modern civilizations evolved, was characterized by relatively 
stable and warm planetary conditions. Human activities have now 
brought Earth outside of the Holocene’s window of environmental 
variability, giving rise to the proposed Anthropocene epoch (4, 5). 
Planetary-scale environmental forcing by humans continues and 
individual Earth system components are, to an increasing extent, in 
disequilibrium in relation to the changing conditions. As a conse- 
quence, the post-Holocene Earth is still evolving, and ultimate 
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global environmental conditions remain uncertain. Paleoclimate re- 
search, however, documents that Earth has previously experienced 
largely ice-free conditions during warm periods (6, 7) with corre- 
spondingly different states of the biosphere. It is clearly in human- 
ity’s interest to avoid perturbing Earth system to a degree that risks 
changing global environmental conditions so markedly. Ice cover is 
only one indicator of substantial system-wide change in numerous 
other Earth system dimensions. The planetary boundaries frame- 
work delineates the biophysical and biochemical systems and pro- 
cesses known to regulate the state of the planet within ranges that 
are historically known and scientifically likely to maintain Earth 
system stability and life-support systems conducive to the human 
welfare and societal development experienced during the Holocene. 

Currently, anthropogenic perturbations of the global environ- 
ment are primarily addressed as if they were separate issues, e.g., 
climate change, biodiversity loss, or pollution. This approach, 
however, ignores these perturbations’ nonlinear interactions and re- 
sulting aggregate effects on the overall state of Earth system. Plane- 
tary boundaries bring a scientific understanding of anthropogenic 
global environmental impacts into a framework that calls for con- 
sidering the state of Earth system as a whole. 

For >3 billion years, interactions between the geosphere (energy 
flow and nonliving materials in Earth and atmosphere) and bio- 
sphere (all living organisms/ecosystems) have controlled global en- 
vironmental conditions. Earth system's state changed in response to 
forcings generated by external perturbations (e.g., solar energy 
input and bolide strikes) or internal processes in the geosphere 
(e.g., plate tectonics and volcanism) or biosphere (e.g., evolution 
of photosynthesis and rise of vascular plants). These forcings were 
processed through interactions and feedbacks among processes and 
systems within Earth system, shaping its often complex overall re- 
sponse. Today, human activities with planetary-scale effects act as 
additional forcing on Earth system. Thus, the anthroposphere has 
become an additional functional component of Earth system (3, 
8), capable of altering Earth system state. The planetary boundaries 
framework formulates limits to the impact of the anthroposphere 
on Earth system by identifying a scientifically based safe operating 
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space for humanity that can safeguard both Earth’s interglacial state 
and its resilience. 

The Holocene state of Earth is the benchmark reference in this 
context, as many of the components comprising the planetary 
boundary framework were rather stable during this period. This is 
also the only Earth system state civilizations have historically 
known. Climate is a manifestation of external forcing, e.g., solar ac- 
tivity, orbital cycles, and interactions among Earth system compo- 
nents, and global mean surface temperature varied by only +0.5°C 
(9) from the Neolithic [~9000 before the present (B.P.)] until the 
Industrial Revolution. Biomes across Earth have also largely been 
stable over the past 10,000 years, with preindustrial global terrestrial 
net primary production (NPP) varying by not >55.9 + 1.1 billion 
tonnes (Gt) of C year! (20) (see the Supplementary Materials). 
Bias-corrected data (10) confirm that preindustrial global precipita- 
tion levels were also stable, particularly from the mid-Holocene 
onward. These data provide strong support for using the Holocene 
(see the Supplementary Materials) as the planetary boundaries ref- 
erence state for a stable and resilient planet. 

All of the framework’'s individual boundaries therefore adopt 
preindustrial Holocene conditions as a reference for assessing the 
magnitude of anthropogenic deviations. Available data and state 
of knowledge from analytics and modeling of the framework com- 
ponents dictate the methods for derivation and quantification of the 
individual boundaries and their precautionary guardrails. Despite 
data constraints, efforts have been made to identify suitable 
control variables for all boundaries, together with evidence of 
how much perturbation leads to generation of impacts or altered 
interactions/feedbacks that can potentially cause irreversible 
changes to Earth's life support systems. The focus is always at 
Earth system rather than regional scale, even when the evidence 
used to establish boundaries originates from regional studies. In 
these cases, regional evidence is combined to assess Earth system 
impacts of cumulative transgressions across multiple region- 
al systems. 

The planetary boundaries framework has attracted considerable 
scientific and societal attention, inspiring governance strategies and 
policies at all levels. The framework evolves through updates made 
in light of recent scientific understanding. Here, we bring together 
advances from different fields of science to update the framework 
and the status of its boundaries. Boundaries are, for the first time, 
proposed for all of the individual components of the framework. 
Updates of the functional biosphere integrity and aerosol loading 
boundaries are based on analyses presented here. Recent analyses 
form the basis for updates of the freshwater change and novel enti- 
ties boundaries. Last, the importance of considering human impacts 
on components of the global environment in a system context is il- 
lustrated using a modeling exercise exploring how various scenarios 
of transgression of the land system (representing the biosphere) and 
climate change boundaries combine to affect Earth system 
characteristics. 


Framework components 

Understanding how biosphere, anthroposphere, and geosphere pro- 
cesses interact with one another is a prerequisite for developing re- 
liable projections of possible future Earth system trajectories. A fully 
process-based understanding of the interactions between these 
domains is, however, still only partially available. The planetary 
boundaries framework calls for more deeply integrated modeling 
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of Earth system by bringing together currently available evidence 
for the relevant processes and their interactions from different dis- 
ciplines and sources. 

The nine boundaries all represent components of Earth system 
critically affected by anthropogenic activities and relevant to Earth's 
overall state. For each of the boundaries, control variables are 
chosen to capture the most important anthropogenic influence at 
the planetary level of the boundary in focus. For example, land 
system change arises from myriad human activities, ultimately ag- 
gregating to alteration of biomes. From a planetary perspective 
however, during the Holocene, forests were the land biome with 
the strongest functional coupling to the climate system (11, 12). 
Therefore, global reduction in forest area is adopted as the control 
variable representing all land system change. Similarly, the control 
variable introduced here for the functional component of the bio- 
sphere integrity boundary, human appropriation of NPP (HANPP), 
focuses on the ability of the biosphere as a whole to provide func- 
tional feedbacks in Earth system. Control variables should ideally 
lend themselves to empirical determination and be computable 
for use in Earth system projections (e.g., process-based simulation 
of future change in forest cover) where possible. 

Boundary positions do not demarcate or predict singular thresh- 
old shifts in Earth system state. They are placed at a level where the 
available evidence suggests that further perturbation of the individ- 
ual process could potentially lead to systemic planetary change by 
altering and fundamentally reshaping the dynamics and spatiotem- 
poral patterns of geosphere-biosphere interactions and their feed- 
backs (13, 14). 

Zone of increasing risk (of Earth system losing Holocene-like 
characteristics) is now used to assess the status for transgressed 
boundaries rather than the “zone of uncertainty" (2) as demarcation 
of this zone is based on more than what is usually referred to as sci- 
entific uncertainty. A large body of recent research [e.g., (15—17)] 
provides strong evidence supporting the conclusion (2) that the 
climate change and biosphere integrity boundaries are in a zone 
of rapidly increasing and systemically linked risks. This strengthens 
the rationale for using the precautionary principle to set the plane- 
tary boundaries at the lower end of the zone of increasing risk. For 
example, for the climate change planetary boundary, we retain the 
boundary of 350 parts per million (ppm) CO, with the zone of in- 
creasing risk ranging from 350 to 450 ppm before reaching high 
risk. This corresponds approximately to a range of global mean 
surface temperature rise of 1° to 2°C (assuming mainstream scenar- 
ios on non-CO, forcing). Precaution places the planetary boundary 
at the start of increasing risk (350 ppm * 1°C), ie., slightly below the 
1.5°C target identified in the Paris Agreement. The 1.5°C target is 
one that science increasingly demonstrates is associated with sub- 
stantial risk of triggering irreversible large change and that crossing 
tipping points cannot be excluded even at lower temperature in- 
creases (18). In recognition of the buffering resilience of Earth 
system, most boundaries are nevertheless set at values higher than 
their observed range through the Holocene up to the Industrial Rev- 
olution (for CO = 280 ppm) (see the Supplementary Materials). 
The stability and characteristic range of variability of interglacial 
Earth system states in Pleistocene paleoclimate (19) and Earth 
system modeling (20) suggest that Earth system would likely 
remain in a stable, Holocene-like state if all boundaries were re- 
spected despite their being at least temporarily outside the envelope 
of Holocene variability. 
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The distinction between zones of “increasing” and “high” risk 
cannot be sharply defined. There is accumulating evidence that 
the current level of boundary transgression has already taken 
Earth system beyond a “safe” zone. However, we still lack a compre- 
hensive, integrated theory, backed by observations and modeling 
studies, that can identify when a transition from a rising level of 
risk to one with very high and dangerous risks of losing a Holo- 
cene-like Earth system state may occur. Therefore, the “burning 
embers” approach introduced by the Intergovernmental Panel on 
Climate Change (IPCC) to represent the gradual transitions from 
moderate (yellow) to high (red) to very high (purple) risks is 
adopted here. 

Throughout Earth's history, geosphere-biosphere interactions 
were an internal driver of Earth system state. The climate change 
planetary boundary is used here as a proxy for the geosphere. There- 
fore, climate change and biosphere integrity are identified as “core 
boundaries” (2) in the framework. The introduction of novel enti- 
ties is a new anthropogenic driver of Earth system change that, if 
sufficiently transgressed, could, on its own, alter Earth system 
state. However, this planetary boundary acts largely through pertur- 
bation of the core boundaries, especially biosphere integrity. In con- 
trast to the definition applied earlier (2) where “naturally occurring 
elements mobilized by anthropogenic activities” were included, the 
definition of novel entities is now restricted to include only entities 
that, in the absence of the anthroposphere, are not present in 
Earth system. 

Quantifying interactions between boundaries remains a major 
challenge. However, some progress has been made since the last 
framework update (2). Recent studies (13, 14, 21, 22) have shown 
that additional or more extensive transgression of one planetary 
boundary can change risk gradients for other boundaries. For 
example, there is increasing evidence to suggest that transgressing 
either the climate change or biosphere integrity planetary boundary 
can potentially lead to more steeply increasing risk in the other (21). 
In the current absence of a comprehensive Earth system model that 
fully captures interactions between all component spheres, we 
explore below how various scenarios of transgression of the land 
system (representing the biosphere) and climate change boundaries 
combine to control biologically mediated carbon storage at the 
planetary level. 


RESULTS 
Biosphere integrity 
Myriad interactions with the geosphere make the biosphere a con- 
stitutional component of Earth system and a major factor in regu- 
lating its state. The planetary functioning of the biosphere 
ultimately rests on its genetic diversity, inherited from natural selec- 
tion not only during its dynamic history of coevolution with the ge- 
osphere but also on its functional role in regulating the state of Earth 
system. Genetic diversity and planetary function, each measured 
through suitable proxies, are therefore the two dimensions that 
form the basis of a planetary boundary for biosphere integrity. As 
applied here, “integrity” does not imply an absence of biosphere 
change but, rather, change that preserves the overall dynamic and 
adaptive character of the biosphere. 

Rockstrém et al. (1) defined the planetary boundary for change 
in genetic diversity as the maximum extinction rate compatible with 
preserving the genetic basis of the biosphere’s ecological 
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complexity. We retain the boundary level of <10 E/MSY (extinc- 
tions per million species-years). The extinction rate control variable 
is challenging to apply in operational contexts, but data and 
methods for directly assessing the genetic diversity component of 
biosphere integrity are emerging [(23) and the Supplementary Ma- 
terials]. Although the baseline rate of extinctions (and of new 
species’ evolution) is both highly variable and difficult to quantify 
with confidence through geological time, the current rate of species 
extinctions is estimated to be at least tens to hundreds of times 
higher than the average rate over the past 10 million years and is 
accelerating (24). We conservatively set the current value for the ex- 
tinction rate at >100 E/MSY (24-26). Of an estimated 8 million 
plant and animal species, around 1 million are threatened with ex- 
tinction (16), and over 10% of genetic diversity of plants and 
animals may have been lost over the past 150 years (23). Thus, the 
genetic component of the biosphere integrity boundary is markedly 
exceeded (Fig. 1 and Table 1). 

Previously, Steffen et al. (2) proposed using the Biodiversity In- 
tactness Index (BII) (27), an empirically based metric of human 
impacts on population abundances, as an interim proxy for func- 
tional biosphere integrity. It was noted, however, that the link of 
BII to Earth system functions remains poorly understood and BII 
cannot be directly linked to the planetary biogeochemical and 
energy flows relevant for establishing Earth system state. In addi- 
tion, BII relies on expert elicitation to estimate temporal changes 
in species abundances/distributions, and this knowledge is not 
readily available for many regions, including the oceans. Martin 
et al. (28) have also recently suggested that BII only partially reflects 
human impacts on Earth system. 

We therefore now replace this metric with a computable proxy 
for photosynthetic energy and materials flow into the biosphere 
(29), ie., net primary production (NPP), and define the functional 
component of the biosphere integrity boundary as a limit to the 
human appropriation of the biosphere’s NPP (HANPP) as a frac- 
tion of its Holocene NPP. NPP is fundamental for both ecosystems 
and human societies as it supports their maintenance, reproduc- 
tion, differentiation, networking, and growth. Biomes depend on 
the energy flow associated with NPP to maintain their planetary 
ecological functions as integral parts of Earth system. NPP-based 
energy flows into human societies should therefore not substantially 
compromise the energy flow to the biosphere (30). The proxy com- 
plements the diversity-based dimensions of biosphere integrity, 
covered by the genetic component, which captures the importance 
of variability in living organisms for the functioning of ecosystems. 
The suitability of NPP and HANPP for defining a planetary boun- 
dary has previously been discussed by Running (31) and Haber] 
et al. (32). 

We determine the terrestrial biosphere’s Holocene NPP to have 
been 55.9 Gt of C year’ (20) and exceedingly stable, varying by not 
more than +1.1 Gt of C year~' despite regional variations in time 
(see the Supplementary Materials). Our model analyses suggest 
that NPP still had a Holocene-like level in 1700 (56.2 Gt of C 
year ' for potential natural vegetation and 54.7 Gt of C year * 
when land use is taken into account). By 2020, potential natural 
NPP would have risen to 71.4 Gt of C year~' because of carbon fer- 
tilization, a disequilibrium response of terrestrial plant physiology 
to anthropogenically increasing CO, concentration in the atmo- 
sphere, whereas actual NPP was 65.8 Gt of C year™' due to the 
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Fig. 1. Current status of control variables for all nine planetary boundaries. Six of the nine boundaries are transgressed. In addition, ocean acidification is approach- 
ing its planetary boundary. The green zone is the safe operating space (below the boundary). Yellow to red represents the zone of increasing risk. Purple indicates the 
high-risk zone where interglacial Earth system conditions are transgressed with high confidence. Values for control variables are normalized so that the origin represents 
mean Holocene conditions and the planetary boundary (lower end of zone of increasing risk, dotted circle) lies at the same radius for all boundaries (except for the 
wedges representing green and blue water, see main text). Wedge lengths are scaled logarithmically. The upper edges of the wedges for the novel entities and the 
genetic diversity component of the biosphere integrity boundaries are blurred either because the upper end of the zone of increasing risk has not yet been quantitatively 
defined (novel entities) or because the current value is known only with great uncertainty (loss of genetic diversity). Both, however, are well outside of the safe operating 


space. Transgression of these boundaries reflects unprecedented human disruption of Earth system but is associated with large scientific uncertainties. 


NPP-reducing effects of global land-use (see the Supplementary 
Materials). 

HANPP designates both the harvesting and the elimination or 
alteration (mostly reduction) of potential natural NPP (32), 
mainly through agriculture, silviculture, and grazing. Terrestrial 
HANPP can be estimated both as a fraction of potential natural 
NPP [15.7% in 1950 and 23.5% in 2020; inferred from (33) and 
the Supplementary Materials] and of Holocene mean NPP (30% 
or 16.8 Gt of C year™' in 2020; see the Supplementary Materials). 
We argue that an NPP-based planetary boundary limiting HANPP 
should be set in relation to preindustrial Holocene mean NPP and 
not the current potential natural NPP. This is because the global in- 
crease in NPP due to anthropogenic carbon fertilization constitutes 
a resilience response of Earth system that dampens the magnitude of 
anthropogenic warming. Hence, the NPP contribution to a carbon 
sink associated with CO, fertilization should be protected and sus- 
tained rather than considered as being available for harvesting. Ex- 
amples of large land areas under human use with declining carbon 
sinks, some even turning into carbon sources, i.e., due to human 
overexploitation of biomass, are already being observed, for 
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example, in some Amazonian regions (34) and northern European 
forests. 

As NPP is the basis for the energy and materials flow that under- 
pins the biosphere’s functioning (30), we argue that today’s plane- 
tary-scale impact of HANPP is reflected in the observation that 
major indicators of the state of the biosphere show large and wor- 
risome declines in recent decades (16). This suggests that current 
HANPP is well beyond a precautionary planetary boundary 
aiming to safeguard the functional integrity of the biosphere and 
likely already into the high-risk zone. We therefore provisionally 
set the functional component of the biosphere integrity planetary 
boundary at human appropriation of 10% of preindustrial Holocene 
mean NPP, shifting into the zone of high risk at 20%. The boundary 
thus defined was transgressed in the late 19th century, a time of con- 
siderable acceleration in land use globally (35) with strong impacts 
on species (27), already leading to early concerns about the effects of 
this large-scale land transformation. 

Thus, while the climate warming problem became evident in the 
1980s, problems arising in functional biosphere integrity due to 
human land use began a century earlier. Since the 1960s, growth 
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Table 1. Current status for the planetary boundaries. 


Earth 
system process 


Control variable(s) 


Planetary boundary 


Preindustrial 
Holocene 
base value 


Upper end of 
zone of 
increasing risk 


Current value of control 
variable 


Climate change 


Atmospheric CO, 
concentration (ppm CO) 


350 ppm CO, 


280 ppm CO, 


450 ppm CO, 


417 ppm CO; (417) 


Total anthropogenic 
radiative forcing at top-of- 
atmosphere (W m7”) 


Change in 
biosphere 
integrity 


<10 E/MSY but with an 
aspirational goal of ca. 1 E/ 
MSY (assumed background 

rate of extinction loss) 


Functional integrity: 
measured as energy 
available to ecosystems 
(NPP) (% HANPP) 


HANPP (in billion tonnes of C 
year” ') <10% of preindustrial 
Holocene NPP, i.e., >90% 
remaining for supporting 
biosphere function 


1.9% (20 
variability of 
preindustrial 

Holocene century- 
mean NPP) 


20% HANPP 


30% HANPP (see the 
Supplementary Materials) 


Stratospheric 
ozone depletion 


Stratospheric O3 
concentration, (global 
average) (DU) 


<5% reduction from 
preindustrial level assessed 
by latitude (~276 DU) 


Ocean 


Carbonate ion 
concentration, average 
global surface ocean 
saturation state with respect 
to aragonite (Oarag) 


280% Oarag Of mean 
preindustrial aragonite 
saturation state of surface 
ocean, including natural diel 
and seasonal variability 


Biogeochemical 
flows: P and 
N cycles 


Phosphate global: P flow 
from freshwater systems into 
the ocean; regional: P flow 
from fertilizers to erodible 
soils (Tg of P year~') 


Phosphate global: 11 Tg of P 
year |; regional: 6.2 Tg of P 
year_' mined and applied to 
erodible (agricultural) soils. 
Boundary is a global average, 
but regional distribution is 
critical for impacts. 


1 


0 Tg of P year™ 


Global: 100 Tg of 
P year"; 
regional: 11.2 Tg 


of P year"! 


Global: 22.6 Tg of P year~' 


(75); regional: 17.5 Tg of P 
year | (76) 


Nitrogen global: industrial 
and intentional fixation of N 
(Tg of N year') 


Nitrogen global: 62 Tg of N 
year~', Boundary is a global 
average. Anthropogenic 
biological N fixation on 
agriculture areas highly 
uncertain but estimates in 
range of ~30 to 70 Tg of N 
year_'. Boundary acts as a 
global “valve” limiting 
introduction of new reactive 
N to Earth system, but 
regional distribution of 
fertilizer N is critical 
for impacts. 


0 Tg of N year! 


190 Tg of N year” (84) 


system change 


Global: area of forested land 
as the percentage of original 
forest cover; biome: area of 
forested land as the 
percentage of potential 
forest (% area remaining) 


Global: 75% values are a 
weighted average of the 
three individual biome 
boundaries; biomes: tropical, 
85%; temperate, 50%; 
boreal: 85% 


Global: 54%; 
biomes: tropical, 
60%; temperate, 

30%; 

boreal: 60% 


Global: 60% [(72, 97) and see 


the Supplementary 

Materials]; tropical: 
Americas, 83.9%; Africa, 

54.3%; Asia, 37.5%; 


temperate: Americas, 51.2%; 
Europe, 34.2%; Asia, 37.9%; 


boreal: Americas, 56.6%; 
Eurasia: 70.3% 


Freshwater 
change 


Blue water: human induced 
disturbance of blue 
water flow 


Upper limit (95th percentile) 
of global land area with 
deviations greater than 

during preindustrial, Blue 
water: 10.2% 


9.4% (median of 
preindustrial 
conditions) 


18.2% (46) 


continued on next page 
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Earth 
system process 


Control variable(s) 


Planetary boundary 


Current value of control 
variable 


Preindustrial 
Holocene 
base value 


Upper end of 
zone of 
increasing risk 


water available to plants (% 
land area with deviations 


9.8% (median of 
preindustrial 


from preindustrial conditions) 
variability) 
Atmospheric Interhemispheric difference 0.1 (mean annual 0.03 0.25 0.076 (55, 57, 68) 
aerosol loading in AOD interhemispheric difference) 
Novel entities Percentage of synthetic 0 0 NA Transgressed 


chemicals released to the 
environment without 
adequate safety testing 


in global population and consumption further accelerated land use, 
driving the system further into the zone of increasing risk. HANPP 
has always sustained humanity's need for food, fiber, and fodder, 
and this will continue to be the case in the future, as well as for sus- 
tainable societies. The NPP required to support future societies 
must, however, increasingly be generated through additional pro- 
duction of NPP above the Holocene baseline, not including the 
NPP generated for biology-based carbon sinks. Feeding 10 billion 
people, for example, is theoretically possible within planetary 
boundaries but requires a number of far-reaching transformations 
to improve the impacts of production and regulate demand (36). 

To develop a deeper foundation for the HANPP-based planetary 
boundary for functional biosphere integrity, we need an improved 
understanding of how ecological dynamics generate the functions of 
the biosphere in Earth system. Analysis of NPP should be spatially 
explicit and augmented by computable metrics of ecological desta- 
bilization due to climate and land use pressures, e.g., a metric of bi- 
ogeochemical disruption (37). 

HANPP can also be quantified for marine systems. About two- 
thirds of the ocean area where HANPP is >10% is found above the 
shallow shelf areas (38) where ecosystems are most intensely ex- 
ploited. Regionally, fish catches exceed thresholds of sustainable ex- 
ploitation (39). However, in contrast to land, where most HANPP 
occurs in the form of plant material, i.e., at the lowest trophic level, 
HANPP in the ocean tends to take place at higher trophic levels. 
This means that while HANPP reduces the absolute amount of 
energy available to higher trophic levels on land, much of the 
energy fixed through NPP is used in marine ecosystems before 
HANPP occurs. When the abundance of organisms at the highest 
trophic levels is reduced, changes in marine ecosystem structure 
may change energy flow in these ecosystems (40). Thus, in the 
marine realm, HANPP likely changes the flows rather than the 
amount of energy available. More information about the impacts 
of HANPP in the marine realm is necessary to integrate consider- 
ation of the marine systems in the functional biosphere integrity 
planetary boundary. 


Climate change 

Climate change control variables and boundary levels are retained 
(1, 2). The most important drivers of anthropogenic impacts on 
Earth's energy budget are the emission of greenhouse gases and 
aerosols, and surface albedo changes (17). The control variables 
in the framework are the annual averages of atmospheric CO, 
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concentration and the change in radiative forcing. The planetary 
boundary for atmospheric CO, concentration is set at 350 ppm 
and for radiative forcing at 1 W m ’. Currently, the estimated 
total anthropogenic effective radiative forcing is 2.91 W m ~~ 
[2022 estimate, relative to 1750 (17)], and atmospheric CO, concen- 
tration is 417 ppm [annual mean marine surface value for 2022 
(41)], i.e., further outside the safe operating space on both measures 
than in the last update (2). The 350-ppm boundary would lead to a 
lower level of anthropogenic global warming than the internation- 
ally agreed 1.5°C target in the United Nations Paris Climate Agree- 
ment but is consistent with recent studies (17, 18, 42) suggesting the 
possibility of extreme Earth system impacts even at 1.5° warming, 
with risks increasing already markedly above 1° warming. 


Novel entities 
The definition of this boundary is now restricted to truly novel an- 
thropogenic introductions to Earth system. These include synthetic 
chemicals and substances (e.g., microplastics, endocrine disruptors, 
and organic pollutants); anthropogenically mobilized radioactive 
materials, including nuclear waste and nuclear weapons; and 
human modification of evolution, genetically modified organisms 
and other direct human interventions in evolutionary processes. 
Novel entities serve as geological markers of the Anthropocene 
(5). However, their impacts on Earth system as a whole remain 
largely unstudied. The planetary boundaries framework is only con- 
cerned with the stability and resilience of Earth system, i.e., not 
human or ecosystem health. Thus, it remains a scientific challenge 
to assess how much loading of novel entities Earth system tolerates 
before irreversibly shifting into a potentially less habitable state. 
Hundreds of thousands of synthetic chemicals are now produced 
and released to the environment. For many substances, the poten- 
tially large and persistent effects on Earth system processes of their 
introduction, particularly on functional biosphere integrity, are not 
well known, and their use is not well regulated. Humanity has re- 
peatedly been surprised by unintended consequences of this release, 
e.g., with respect to the release of insecticides such as DDT and the 
effect of chlorofluorocarbons (CFCs) on the ozone layer. For this 
class of novel entities, then, the only truly safe operating space 
that can ensure maintained Holocene-like conditions is one 
where these entities are absent unless their potential impacts with 
respect to Earth system have been thoroughly evaluated. This would 
imply that the quantified planetary boundary should be set at zero 
release of synthetic chemical compounds to the open environment 
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unless they have been certified as harmless and are monitored. That 
is the target set by the Montreal Protocol with respect to the sub- 
stances shown to be harmful by contributing to depletion of the 
ozone layer. 

In their analysis of various strategies for establishing a planetary 
boundary for novel entities, Persson et al. (43) identified the share 
of released chemicals with adequate safety assessment and monitor- 
ing as a candidate control variable. We here adopt this metric. The 
planetary boundary is then set at the release into Earth system of 0% 
of untested synthetics. When synthetics released to the environment 
are thoroughly tested, the ensuing risk of damaging effects is 
lowered. Admittedly, this approach has weaknesses: Data availabil- 
ity is incomplete; safety studies often focus on narrowly defined tox- 
icity and do not capture the “cocktail effects" of chemical mixtures 
in the environment nor their effects under specific conditions. The 
percentage of untested synthetics released globally is unknown. 
However, Persson et al. (43) report that for the chemicals currently 
registered under the EU Registration, Evaluation, Authorisation and 
Restriction of Chemicals (REACH) regulation (a small subset of the 
chemical universe), ~80% of these chemicals had been in use for at 
least 10 years without yet having undergone a safety assessment. 
Likewise, few safety studies consider potential Earth system 
effects. With such an enormous percentage of untested chemicals 
being released to the environment, a novel entities boundary 
defined in this manner is clearly breached. Persson et al. (43) did 
not identify or quantify a singular planetary boundary for novel en- 
tities but, nevertheless, also concluded that the safe operating space 
is currently overstepped. 


Stratospheric ozone depletion 

Stratospheric ozone depletion is a special case related to the anthro- 
pogenic release of novel entities where gaseous halocarbon com- 
pounds from industry and other human activities released into 
the atmosphere lead to long-lasting depletion of Earth's ozone 
layer. The boundary for the safe operating space is set at 276 
Dobson units (DU), i-e., allowing a <5% reduction from the prein- 
dustrial level of 290 DU, assessed by latitude (1). Following the rat- 
ification of the Montreal Protocol in 1987, the trend and global 
extent of ozone depletion have recovered slightly (44, 45). The 
current (2020) global estimate is 284 DU (see the Supplementary 
Materials). Thus, the human perturbation of the stratospheric 
ozone depletion has decreased and is now within the safe operating 
space. The boundary for ozone depletion is currently only trans- 
gressed over the Antarctic and southern high latitudes and only 
in the 3-month Austral spring (45). 


Freshwater change 

To comprehensively reflect anthropogenic modifications of Earth 
system functions of freshwater, this boundary is revised to consider 
changes across the entire water cycle over land (46-48). We here use 
streamflow as a proxy to represent blue water (surface and ground- 
water) and root-zone soil moisture to represent green water (plant- 
available water) (46-48). Control variables are defined as the per- 
centage of annual global ice-free land area with streamflow/root- 
zone soil moisture deviations from preindustrial variability (46, 
48). The new green water component directly accounts for hydro- 
logical regulation of terrestrial ecosystems, climate, and biogeo- 
chemical processes (48), whereas the blue water component 
accounts for river regulation and aquatic ecosystem integrity (46). 
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Moreover, this boundary now captures Earth system impacts of 
both water increases and decreases on a monthly scale and includes 
their spatial patterns (see the Supplementary Materials). 

The control variables describe deviations from the preindustrial 
(here, 1661-1860) state, first determined at the 30 arc-min grid cell 
scale and further aggregated to a global annual value. For both blue 
and green water control variables, boundaries are set at the 95th per- 
centile of preindustrial variability, i.e., variability of the percentage 
of global land area with deviations [~10% for blue and ~11% for 
green water; (46) and the Supplementary Materials]. We assume 
that preindustrial conditions are representative of longer-term Ho- 
locene conditions and that notable deviation from this state puts 
freshwater's Earth system functions at risk. Pending comprehensive 
assessment of impacts of different transgression levels of the blue 
and green water boundaries (e.g., reduced carbon sequestration ca- 
pacity, climate regulation, and biodiversity loss; see the Supplemen- 
tary Materials), the boundary settings are preliminary and highly 
precautionary. Currently, ~18% (blue water) and ~16% (green 
water) of the global land area experience wet or dry freshwater de- 
viations (46). Thus, in contrast to the earlier planetary boundary as- 
sessments (1, 2) where only blue water removal was considered, this 
new approach indicates substantial transgression of the freshwater 
change boundary. Transgressions of both the blue and green water 
boundaries occurred a century ago, in 1905 and 1929, respectively 
(46). Thus, with the revised definition of the control variables, fresh 
water would have been considered transgressed already at the time 
of the previous planetary boundary assessments. The previous 
global-scale control variable would still indicate freshwater use to 
remain in the safe zone, even with newer data sources than those 
used in (1, 2). Recent estimates of global blue water consumption 
totals ~1700 km? year”! (49), i.e., far below the previous boundary 


set at 4000 km? year™’. 


Atmospheric aerosol loading 

Aerosols have multiple physical, biogeochemical, and biological 
effects in Earth system, motivating their inclusion as a planetary 
boundary (see the Supplementary Materials). Anthropogenic 
aerosol loading has increased (50). Changes since the preindustrial 
for natural aerosols (e.g., desert dust, soot from wildfires) are diffi- 
cult to assess because of model differences in the sign of trends (51), 
but observational evidence suggests a global doubling of dust dep- 
osition since 1750 (52). At present, the Sahara is the world's largest 
dust source region [e.g., (53)], but earlier in the Holocene, it was a 
vegetated landscape with many lakes and wetlands (14,500 to 5000 
B.P.). Changes in monsoon rainfalls, involving vegetation-dust- 
climate feedbacks, are thought to have terminated the “green 
Sahara,” leading to major displacements of human settlements 
across parts of Africa and Asia (54). 

Quantification of the aerosol loading planetary boundary is 
hampered by their multiple natural and human-caused sources, dif- 
ferences in chemical composition, seasonality and atmospheric life- 
times, and the consequently very large spatial and temporal 
heterogeneity in distribution and climatic and ecological impacts 
of aerosols. Nevertheless, aerosol optical depth (AOD) provides a 
generic control variable for aerosol loading. AOD is an integrated 
measure of the overall reduction in sunlight reaching Earth's 
surface caused by all absorption and scattering in the vertical air 
column. On the basis of the evidence of the impacts of large AOD 
on regional precipitation over southern Asia, Steffen et al. (2) set a 


7 of 16 


€ZOT ‘ET Joquiajdag uo B1o‘a0uaTOS’ MM M//:Sd}IY WOIT papeo[UMog 


SCIENCE ADVANCES | RESEARCH ARTICLE 


provisional regional planetary boundary of AOD = 0.25 (0.25 to 0.5) 
on the basis that higher AOD values in monsoon regions likely lead 
to significantly lower rainfall, ultimately affecting biosphere integ- 
rity. The annual mean AOD in southern Asia is currently about 0.3 
to 0.35 (55, 56). The current value for the East China region is 0.4 
(55). Thus, aerosol loading in these regions has likely exceeded the 
regionally defined boundary, but with high uncertainty. Data and 
assessments of aerosol impacts on climate and ecosystems are 
lacking to determine whether this regionally defined boundary is 
applicable elsewhere. Global mean AOD at present is 0.14 (57), 
with much higher levels in some regions and with very strong gra- 
dients from land to open ocean (56). 

In addition to the direct effects of AOD on regional climate and 
precipitation, asymmetries in AOD between northern and southern 
hemispheres can affect multiple monsoon systems, as seen for the 
West African monsoon (58) and Indian monsoon (59, 60). The in- 
terhemispheric difference in AOD affects regional monsoon rainfall 
by shifting the location of the Intertropical Convergence Zone (61). 
Large asymmetries in the temperature of northern and southern 
hemispheres arise from differences in natural and anthropogenic 
aerosol emissions, land cover, and other climate forcers (58, 59, 
62, 63). The asymmetric radiative forcing resulting from aerosol 
effects leads to a relative cooling of the northern hemisphere and 
a southward shift in tropical precipitation (64). The interhemi- 
spheric AOD difference and its impact on tropical precipitation 
and water availability are sensitive to the particle size and latitudinal 
and altitudinal distribution of aerosols (65). Studies of aerosol- 
climate interactions following volcanic eruptions (66) indicate 
that monsoon precipitation in the northern hemisphere is weak- 
ened when northern hemisphere AOD is higher and the interhemi- 
spheric AOD difference is greater and is enhanced when more 
aerosols are emitted in the southern hemisphere (smaller inter- 
hemispheric AOD difference). This understanding is broadly con- 
sistent with the decrease in tropical mean precipitation after major 
volcanic eruptions in observations and global climate models (67). 
The IPCC ARé has assessed that observed decreases in global land 
monsoon precipitation from the 1950s to the 1980s are partly attrib- 
uted to human-caused northern hemisphere aerosol emissions, thus 
relatively larger interhemispheric difference (17). In addition to vol- 
canic aerosols, monsoon dynamics and the associated regional rain- 
falls also respond to changes in anthropogenic aerosols (see the 
Supplementary Materials). 

We therefore propose the annual mean interhemispheric differ- 
ence in AOD as a globally defined control variable for aerosol 
loading. The present-day interhemispheric difference is ~0.076 + 
0.006 (mean + SD), based on 12 observational estimates, reaching 
~0.1 in the boreal spring and summers, due to the seasonal increase 
in dust storms that dominate in the northern hemisphere (55). The 
preindustrial annual mean value is estimated as ~0.03, based on 
multimodel analyses (68), indicating an increase in interhemispher- 
ic AOD difference by ~0.04 in the industrial era. Present-day inter- 
hemispheric AOD difference is consistent with Coupled Model 
Intercomparison Project 6 (CMIP6) emission inventories that 
show more anthropogenic aerosols in the northern hemisphere, 
with future projections suggesting a decrease in the asymmetry (69). 

We assign a planetary boundary value of 0.1 for the mean annual 
interhemispheric difference in AOD, with high uncertainty about 
the zone of increasing risks, 0.1 to 0.25. In setting this boundary, 
we note that the impacts of aerosol loading on tropical monsoon 
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systems are already seen today, and the impact is not only restricted 
to rainfall but also affects regional climate more broadly. Aerosol- 
cloud interaction might exacerbate effects of AOD asymmetry. The 
contribution of aerosol-cloud interactions to the hemispheric asym- 
metry of reflected shortwave radiation is unclear. Take for instance 
the current range of anthropogenic aerosol effective radiative 
forcing for present day that has been reported to be —1.6 to —0.6 
W m ” in the global mean for the 16 to 84% confidence interval, 
with aerosol-cloud interactions as a major source for uncertainty 
(51). Other large-scale effects of aerosols, such as air quality 
impacts on land and marine ecosystems, are also already evident 
(17, 70). Biogenic aerosols have not been considered, despite their 
role in feedbacks in Earth system. A much better systemic and quan- 
titative understanding of the hydroclimatic, ecological, and biogeo- 
chemical effects of asymmetric aerosol forcing is needed to refine 
the aerosol loading boundary. 


Ocean acidification 

The control variable used is the carbonate ion concentration in 
surface seawater (specifically, OQarag, the average global surface 
ocean saturation state with respect to aragonite). The original boun- 
dary quantification [>80% of the preindustrial averaged global O.arag 
of 3.44 (1)] is retained. A recent estimate sets the current Olarag at 
~2.8 (71) (see the Supplementary Materials), approximately 81% 
of the preindustrial value. Thus, anthropogenic ocean acidification 
currently lies at the margin of the safe operating space, and the trend 
is worsening as anthropogenic CO, emission continues to rise. 


Land system change 

This boundary focuses on the three major forest biomes that 
globally play the largest role in driving biogeophysical processes 
(2), i.e. tropical, temperate, and boreal. The control variable 
remains the same: forest cover remaining compared to the potential 
area of forest in the Holocene (2). The boundary positions remain at 
85%/50%/85% for boreal/temperate/tropical forests (cf. Table 1 and 
the Supplementary Materials). On the basis of 2019 land-cover 
classification maps derived from satellite observations (72), the 
current state of the regional biomes is similar to that in 2015 
although, for most regions, the amount of deforestation has in- 
creased since 2015 (see the Supplementary Materials). Land-use 
conversion and fires are causing rapid change in forest area (73, 
74), and deforestation of the Amazon tropical forest has increased 
such that it has now transgressed the planetary boundary (Table 1). 
Changes in the methodology and technology used to estimate forest 
cover since 2015 may be influencing the biome-level differences re- 
ported here compared to the last update (2). Nevertheless, there is 
little doubt that the global forest area continues to decrease (74). 


Biogeochemical flows 

Biogeochemical flows reflect anthropogenic perturbation of global 
element cycles. Currently, the framework considers nitrogen (N) 
and phosphorus (P) as these two elements constitute fundamental 
building blocks of life, and their global cycles have been markedly 
altered through agriculture and industry. Anthropogenic impacts 
on global carbon cycling are equally fundamental but are addressed 
in the climate and biosphere integrity boundaries. Other elements 
could come into focus under this boundary as an understanding of 
human perturbation of element cycles advances. For both N and P, 
the anthropogenic release of reactive forms to land and oceans is of 
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interest, as altered nutrient flows and element ratios have profound 
effects on ecosystem composition and long-term Earth system 
effects. Some of today's changes will only be seen on evolutionary 
time scales, while others are already affecting climate and biosphere 
integrity. 

For P, we retain the regional-level and global boundaries pro- 
posed by Steffen et al. (2). The global boundary for P is a sustained 
flow of 11 Tg of P year~’ from fresh water to the ocean, to avoid 
large-scale anoxia. We have not found newer studies quantifying 
P flows in fresh water to the sea since that used for the 2015 frame- 
work update, ie., an estimated 22 Tg of P year ' (75). The regional 
level boundary is set at a flow of 6.2 Tg of P year™' from fertilizers to 
erodible soils, to avert widespread eutrophication of freshwater eco- 
systems. The current rate of application of P in fertilizers to crop- 
lands is 17.5 Tg of P year‘ (76) although P use is rising and much 
higher estimates of up to 32.5 Tg of P year * have been reported in 
other studies (77-79). Thus, both the global and regional boundar- 
ies for P are exceeded. The planetary boundary for N is the applica- 
tion rate of intentionally fixed N to the agricultural system of 62 Tg 
of N year~' [unchanged from (2)]. Currently, the application of in- 
dustrially fixed N fertilizer is 112 Tg of N year! (80). Quantifica- 
tion of anthropogenic biological N fixation in connection with 
agriculture is highly uncertain, but the most recent estimates are 
in the range of ~30 to 70 Tg of N year! (81-83). According to 
Food and Agriculture Organization (84), the total introduction of 
anthropogenically fixed N applied to the agricultural system is 
~190 Tg year‘ so this boundary is also globally transgressed. 


DISCUSSION 

Six planetary boundaries are found currently to be transgressed 
(Fig. 1 and Table 1). For all of the boundaries previously identified 
as transgressed [climate change, biosphere integrity (genetic 
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diversity), land system change, and biogeochemical flows (N and 
P)], the degree of transgression has increased since 2015. We have 
introduced HANPP as a control variable for the functional compo- 
nent of biosphere integrity and argue that this boundary is also 
transgressed. Drawing on the considerable recent scientific progress 
made in refining the safe operating space for water, control variables 
for both green and blue water components are now included in the 
freshwater change planetary boundary. The boundary is trans- 
gressed for both components. Global boundaries for aerosol 
loading and novel entities are proposed. The novel entities boun- 
dary is transgressed. The global aerosol loading boundary is not 
transgressed although regional transgressions are noted. 


Earth system effects of differing scenarios of transgression 
of land system change and climate boundaries 

To illustrate the importance of considering the multiple anthropo- 
genic impacts on the global environment in a systemic context 
rather than individually, we examine how varying degrees of trans- 
gression of the climate and land system change boundaries combine 
to influence two codeterminants of Earth system state: temperature 
and terrestrial carbon storage. 

For climate change, the Potsdam Earth Model (POEM) [(85) and 
the Supplementary Materials] is forced by increased atmospheric 
carbon dioxide levels (350, 450, and 550 ppm), and land system 
change is forced with land-use patterns representing different 
extents of tropical, temperate, and boreal forest cover (see the Sup- 
plementary Materials). As some biological processes take centuries 
to approach a steady state, we investigate changes in both the short 
(1988-2100) and the long term (2100-2770). This also enables us to 
examine the veracity of the placement of these planetary boundaries 
and their zones of increasing risk in terms of critical Earth system 
responses. 


@ Beyond zone 
of increasing 
risk: 40/20/40% 


@ On upper end of 
zone of increasing 
risk: 60/30/60% 


OShort @ Long 
term term 
(2100) @ (2770) 


+200 Gt of C 


Land C Source 
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Fig. 2. Impact of the combined effect of land system change and climate change boundary states on trajectories of terrestrial carbon stocks and global land 
temperature. Results are based on idealized Earth system model experiments with varying planetary boundary status, ranging from maintaining the planetary boundary 
(85%/50%/85% boreal/temperate/tropical forest remaining, 350-ppm atmospheric CO3, green), the upper end of the zone of increasing risk (60%/30%/60%, 450 ppm, 
orange), and beyond the zone of increasing risk (40%/20%/40%, 550 ppm, red). Open circles represent the short-term changes (1988-2100) of the system, while colored 
circles the long-term changes (2100-2770). Their colors denote the state of the land system change boundary, while the climate change boundary is shown on the y axis. 
The locations of the circles on the x axis represent the changes in the land carbon stocks, and the associated land temperature changes are given next to each circle, both 
compared to the year 1988. Transgressing the climate change boundary (y axis) is mostly connected to an increase in temperature, while the transgression of land system 
change leads to a loss of terrestrial carbon stocks (source) of 100 to 200 Gt of C. 
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According to these simulations, anthropogenic activities 
brought both climate and land system change outside of their safe 
operating space around 1988. Had Earth system remained forced by 
1988 conditions (350 ppm and 85%/50%/85% of tropical/temper- 
ate/boreal forest cover remaining), the simulations show that tem- 
perature over the global land surface would not have increased by 
more than an additional 0.6°C in the subsequent 800 years (and not 
>1.3°C compared to the preindustrial period). Only a small (cumu- 
lative 25 Gt of C) terrestrial carbon source would have developed by 
2100 and a cumulative source of not >68 Gt of C after 800 years. 
Thus, the exercise suggests that essentially stable planetary condi- 
tions would have been maintained had human impacts on these 
two boundaries remained at their 1988 levels, i.e., marginally 
within the safe operating space. 

Both of these planetary boundaries have, however, since been 
transgressed into a zone of increasing risk of systemic disruption. 
If climate and land system change can be halted at 450 ppm and 
forest cover retained at 60%/30%/60% of boreal/temperate/tropical 
natural cover, then the simulation indicates a mean temperature rise 
over land of 1.4°C by 2100 (in addition to 0.7°C between preindus- 
trial time and 1988) and 1.9°C after 800 years as vegetation evolves 
in a warmer climate and associated carbon fertilization (Fig. 2). 

Carbon fertilization of vegetation growth counters the negative 
impacts of climate warming on the global average carbon sinks, 
leading to only moderate cumulative loss in terrestrial carbon due 
to additional deforestation. If, however, deforestation had been 
maintained at the level of the planetary boundary rather than 
having been allowed to rise in the zone of increasing risk, then 
the land biosphere would have developed a cumulative carbon 
sink rather than a source, contributing to stabilizing Earth's condi- 
tions. In contrast, if deforestation is allowed to breach into the high- 
risk zone, then simulations show a substantial additional carbon 
leakage to the atmosphere both over the short and long term (132 
and 211 Pg of C), despite strong CO, fertilization of vegetation 
growth in the model (Fig. 2). 

The situation is even more extreme if atmospheric CO, concen- 
tration rises above the risk zone (550 ppm; Fig. 2) and deforestation 
continues. Not only is the temperature on land about 2.7°C warmer 
than in 1988 (3.4°C warmer than preindustrial), but also around 145 
Gt of C would be lost long-term from terrestrial vegetation and soils. 
Note that these findings reflect optimistic modeling assumptions on 
carbon fertilization. Many of the ecological factors not sufficiently 
represented in current biogeochemical models could lead to even 
less desirable consequences of leaving the safe operating space. 
These simulations illustrate clearly that human impacts on 
climate and forest cover must be considered in a systemic context. 


They furthermore support the placement of the planetary boundar- 
ies for climate and land system change at the lower end of the zone 
of increasing risk. 


Influence of climate change on biologically mediated C 
sinks in the ocean 

Approximately 450 Gt of C is bound up in terrestrial biota, primar- 
ily in plants (86), while only ~6 Gt of C is found in ocean biota (87). 
Biologically mediated marine carbon sinks are composed of partic- 
ulate organic carbon (POC) that can potentially sink below the per- 
manent thermocline (biological pump) and dissolved organic 
C. Via microbial breakdown of POC and dissolved organic C, 
CO, is released. When this release influences partial pressure of 
CO, in surface waters, it tends to reduce oceanic carbon uptake 
from the atmosphere. Microbial respiration is highly sensitive to 
temperature and, in a warmer ocean, an increased release of CO, 
in surface waters is predicted (88). The biologically mediated 
carbon sink in the ocean most exposed to climate change is the 
amount of carbon fixed by photosynthesis (NPP), i.e., POC, in 
the surface ocean that is ultimately transported into the ocean inte- 
rior via the biological pump. When this occurs, the resulting carbon 
drawdown reduces partial pressure of CO, in the surface layer and 
tends to increase the atmosphere-to-ocean CO) flux. 

These biological processes are implicitly and, in some cases, ex- 
plicitly included in the CMIP6 models informing the IPCC. 
However, as these models configure biologically mediated carbon 
flows differently, there is considerable variability in their results. 
Models used by the IPCC do not even agree on the direction of 
change in NPP in response to climate change (89). Our model 
runs (see the Supplementary Materials) suggest no significant 
change in globally averaged ocean NPP under the different 
climate forcing conditions and only a modest decrease in exported 
material out of the surface layer [new production (ANP); Table 2]. 
Using empirical relationships (90, 91) describing the transfer of 
carbon to the ocean interior and derived from the contemporary 
ocean to estimate biological pump sensitivity to future temperature 
increases indicates a similar weakening of the pump in the upper 
ocean (Table 2 and the Supplementary Materials). That these two 
independent methods indicate similar decreases in the export of 
POC from the surface layer lends confidence both in the direction 
and magnitude of climate impacts on this biologically mediated 
global carbon sink. 

The analysis shows that DIC (dissolved inorganic carbon; in- 
cluding CO,) accumulates over time in the ocean as a whole, par- 
ticularly in the upper ocean (<1000 m; Table 2). Changes in the 
biologically driven accumulation rates are relatively small compared 


Table 2. Global averaged change in three scenarios from the initial state (1988-2018): change in sea surface temperature (ASST), new production (ANP), 
and biogenic particulate flux below 500 m depth (AFso0,,) including model and empirically derived values, surface saturation state of aragonite (AQ), 
and the DIC inventory between the surface and 1000 m depth (ADIC,_jo00m)- 


Scenario ASST ANP AFs500m AQ ADICo_1000m 

(ppm) (°C) Model (%) Empirical (%) Model (%) Empirical (%) (-) (Gt of C) 

350 0.3 2.0 2.5 1.9 1.8 0.0 38 

450 1.0 0.0 14 0.0 -3.5 -0.4 172 

550 e/ —2.5 —1.0 -3.1 -9.4 —0.7 273 
Richardson et al., Sci. Adv. 9, eadh2458 (2023) 13 September 2023 10 of 16 


€ZOT ‘ET Joquiajdag uo B10'a0UNTOS’ MM M//:Sd}IY WOIT papeo[UMog 


SCIENCE ADVANCES | RESEARCH ARTICLE 


to the change in the total DIC inventory that is mainly driven by the 
solubility pump, i-e., the tendency of increased oceanic uptake when 
atmospheric partial pressure of CO) rises. The organic matter flux 
below the 500-m-depth horizon (AFso9) varies between 3 and 9% 
between the model and empirically derived fluxes in the 550 ppm 
scenario with the model-derived sensitivity being lowest. This illus- 
trates the current uncertainty in quantifying climate-driven feed- 
backs on the biological pump. The implied accumulation of DIC 
in the surface ocean will tend to decrease the uptake of atmospheric 
CO,, thus counteracting global actions for stabilizing or even reduc- 
ing atmospheric CO, concentrations. The ocean response to 
reduced greenhouse gases will be complex and occur on different 
time scales, e.g., the characteristic response time simulated for the 
total carbon pool in the upper 1000 m is ~150 years (550 ppm; see 
the Supplementary Materials). However, the natural ocean carbon 
sink will gradually decrease on millennial time scales. 

The reduction of sinking organic material will affect the meso- 
pelagic ecosystem (i.e., the subsurface ecosystem between 200 and 
1000 m in depth, one of the largest biomes on Earth and one that 
hosts numerous transient grazers, including some whales). The flux 
of organic material via sinking represents the energy source for or- 
ganisms in this biome. A reduction of up to ~10% of energy flux 
would potentially have enormous consequences for this biome 
and, thereby, its biosphere integrity. Recent paleontological recon- 
structions (92) provide evidence that these decreases in carbon flux 
to the mesopelagic may have occurred in relation to past 
climate changes. 

Acidification due to increased CO; reduces the saturation state of 
aragonite (Q). It tends to hinder the biological formation of calcium 
carbonate, an essential component for shell and reef-forming or- 
ganisms. The relatively short equilibration time of the surface 
ocean with atmospheric CO, implies a response time of © to in- 
creased CO, of only a few decades, comparable to the current acid- 
ification rate (see the Supplementary Materials). The current rate is 
probably a hundred times faster than at any time during the last 
hundreds of millennia (93), confirming the tied relations to trans- 
gression of the climate change boundary, leading to the rising risk of 
weakening ocean biosphere integrity, and worsening the aragonite 
saturation state of the ocean acidification boundary. 


A systemic framework for addressing global anthropogenic 
impacts on Earth system 

The scientific updates and analyses presented here confirm that hu- 
manity is today placing unprecedented pressure on Earth system. 
Perhaps most worrying in terms of maintaining Earth system in a 
Holocene-like interglacial state is that all the biosphere-related plan- 
etary boundary processes providing the resilience (capacity to 
dampen disturbance) of Earth system are at or close to a high-risk 
level of transgression. In a recent study (18), it was shown that 
several regional climate tipping points, relevant for stabilizing the 
global system, have already been or are close to being transgressed, 
thus weakening global resilience capacity. This implies low/falling 
resilience precisely when planetary resilience is needed more than 
ever to cope with increasing anthropogenic disturbances. There is 
an urgent need for more powerful scientific and policy tools for an- 
alyzing the whole of the integrated Earth system with reliability and 
regularity and guiding political processes to prevent altering the 
state of Earth system beyond levels tolerable for today's societies. 
In addition to more consistent collection and collation of relevant 
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global environmental data, this will require the development of 
Earth system models that more completely capture geosphere-bio- 
sphere-anthroposphere interactions than is the case today. The 
known interdependence of planetary boundaries is confirmed by 
Earth system science understanding (14, 22) of the planet as an in- 
tegrated, partially self-regulating, system. To better understand the 
risk to this system and the critical boundaries that humankind 
should consider in its economic and social activities, Earth system 
analysis now has to continue advancing a planetary boundaries 
framework. In addition, it must substantially increase the ecological 
realism of simulation and analyses of the biosphere as an adaptive 
core entity of Earth system. These initiatives are underway but have 
to be further developed into a coherent process of integrated Earth 
system analysis across the physical, chemical, and biological 
domains not focused just on climate. 

Successfully addressing anthropogenic climate change will 
require consideration of internal biosphere-geosphere interactions 
within Earth system. Our model results demonstrate that one of the 
most powerful means that humanity has at its disposal to combat 
climate change is respecting the land system change boundary. 
Bringing total global forest cover back to the levels of the late 
20th century would provide a substantial cumulative sink for atmo- 
spheric CO, in 2100. This reforestation seems unlikely, however, 
given the current focus on biomass as a replacement for fossil 
fuels and the creation of negative CO, emissions via bioenergy 
with carbon capture and storage. Both activities are already 
serving to increase pressure on Earth's remaining forest area. Nev- 
ertheless, our study indicates that failure to respect the land system 
change planetary boundary can potentially jeopardize efforts to 
achieve the global climate goals adopted in the Paris Agreement. 

Meanwhile, this update of the planetary boundaries framework 
may serve as a renewed wake-up call to humankind that Earth is in 
danger of leaving its Holocene-like state. It may also contribute to 
guiding the substantial human opportunities for sustainable devel- 
opment on our planet. Scientific insight into planetary boundaries 
does not limit, but stimulates, humankind to innovation toward a 
future in which Earth system stability is fundamentally preserved 
and safeguarded. 


MATERIALS AND METHODS 
To quantify the aerosol boundary, we consider cases where a natural 
pulse of sulfate aerosol emissions from volcanic eruptions in the 
northern hemisphere led to subsequent rainfall deficits in the 
Sahel. The eruption of El Chichon led to a peak interhemispheric 
AOD difference of 0.07 and that of Katmai to an AOD difference 
of 0.08 (55). We also consider a model study of intentional sulfate 
injections into the stratosphere. This study is based on stratospheric 
aerosols, which have no direct interaction with clouds and vegeta- 
tion. However, it does indicate that an interhemispheric sulfate 
AOD difference of ~0.2 would decrease tropical monsoon precipi- 
tation in the northern hemisphere by ~10% and India's mean pre- 
cipitation by >20% (59). Together, these studies suggest that a raised 
interhemispheric AOD difference caused by persistent and widely 
distributed aerosol emissions could lead to major reductions in pre- 
cipitation in the tropics. 

To examine differing scenarios of transgression of land system 
and climate change boundaries, we use the POEM [(85) and the 
Supplementary Materials], which links models of atmospheric 
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and ocean circulation with models of the marine (BLING) (94) and 
terrestrial biosphere (LPJmL5) [(95) and the Supplementary Mate- 
rials]. We study scenarios where each of these two planetary boun- 
dary dimensions are either fixed at the value of the boundary, a 
value in the zone of increasing risk, or a value in the high-risk 
zone. Once the respective scenario condition is attained, the associ- 
ated level of scenario forcing remains constant, while the long-term 
implications under these fixed conditions evolve. Correspondingly, 
vegetation dynamics (e.g., biome distributions) and related carbon 
pools and fluxes develop according to biophysical climate interac- 
tions under the given forcing conditions, while biogeochemical 
feedbacks on the atmosphere are not considered because of the re- 
spective boundary or transgression forcing remaining fixed. 


Supplementary Materials 
This PDF file includes: 
Supplementary Information 

Figs. S1 to S9 

Tables $1 to $3 
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